Ground-state energy and depletions for a dilute binary Bose gas 
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When calculating the ground-state energy of a weakly interacting Bose gas with the help of the 
customary contact pseudopotential, one meets an artifical ultraviolet divergence which is caused by 
the incorrect treatment of the true interparticle interactions at small distances. We argue that this 
problem can be avoided by retaining the actual, momentum-dependent interaction matrix elements, 
and use this insight for computing both the ground-state energy and the depletions of a binary Bose 
gas mixture. Even when considering the experimentally relevant case of equal masses of both species, 
the resulting expressions are quite involved, and no straightforward generalizations of the known 
single-species formulas. On the other hand, we demonstrate in detail how these latter formulas are 
recovered from our two-species results in the limit of vanishing interspecies interaction. 
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I. INTRODUCTION 

Recently, the experimental investigation of multicom- 
ponent Bose-Einstein condensates has made substantial 
progress P, S S S 13 ■ A considerable amount of the- 
oretical work has been devoted to such systems, focus- 
ing, among others, on mean-field descriptions of trapped 
binary mixtures |a, H, stability properties and phase 
segregation [1^ lialll El El 111, c ollective excita- 
tions UJ,UJ,UJ, spinor condensates and multi- 
component condensates in optical lattices 20,, JlJ . 

In this paper we determine the ground-state energy 
and the depletions for a homogeneous binary conden- 
sate within the Bogoliubov approximation . The cal- 
culation of the ground-state energy beyond the mean- 
field approximation, and of the depletion, of a single- 
species condensate is a major achievement of the the- 
ory of weakly interacting Bose gases designed in the late 
fifties m HI i25, ^ Recent significant devel- 

opments include the first mathematically rigorous proof 
of the leading term of the commonly accepted expression 
for the ground-state energy per particle ,29] , the explicit 
calculation of nonuniversal contributions due to three- 
body scattering within the framework of effective field 
theory m EH N, and the clarification of the relation be- 
tween the Gross-Pitaevskii and Bogoliubov descriptions 
of a dilute Bose gas [S^ • One of the technical difficulties 
presented by the subject stems from the fact that typical 
atomic interaction potentials exhibit a strong increase at 
short distances, so that their Fourier transforms diverge 
for vanishing momentum transfer. Hence, one customar- 
ily replaces the true interaction potential by the contact 
pseudopotential 
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where m denotes the atomic mass, and a is the s- 
wave scattering length for collisions among the cold 
atoms. Thus, one relies on the knowledge that the low- 
temperature properties of dilute Bose gases are governed 
solely by the scattering lengths, an assertion commonly 
referred to as the "Landau postulate" Indeed, the 

leading terms of an expansion of the energy density for 
the ground state of a dilute homogeneous single-species 
gas of Bose particles interacting via some potential with 
scattering length a can be written as |32| 
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first derived by Lee and Yang for the particular case of 
hard-sphere Bosons with diameter a |23L l2l| . Nonuni- 
versal, i. e. potential-dependent contributions figure only 
as higher-order corrections, due to the effective rarige for 
s-wave scattering, and three-body interaction 

While the standard contact pseudopotential ^ is suf- 
ficient for determining the quasiparticle energy disper- 
sion for a weakly interacting Bose gas, problems arise in 
the calculation of its ground-state energy. Namely, when 
coupling bare-particles states with opposite momenta hk 
and — hk for setting up the harmonic oscillator-like states 
corresponding to the phonon quasiparticles, one also has 
to account for the "vacuum" energies of these oscillators. 
The energy of the ground state of the weakly interacting 
gas, defined by the absence of phonons for all k, then ac- 
quires a contribution from the sum over all these vacuum 
energies. However, if one naively uses the pseudopoten- 
tial (Q, all interaction matrix elements are replaced by 
a constant, instead of being properly suppressed at high 
momentum transfer, so that this sum becomes ultraviolet 
divergent. It is understood that this divergence is not a 
fundamental physical problem, but rather a consequence 
of the incorrect handling of the true interparticle interac- 
tion at short distances through the delta- function pseu- 
dopotential : Actually, the proper pseudopotential op- 
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erator derived by Huang and Yang |34| is given by 



II. THE HAMILTONIAN IN QUADRATIC 
APPROXIMATION 
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when operating on wave functions with a 1/r-type sin- 
gularity, this operator yields results quite different from 
those given by its simplification . The modern solution 
to the problem is based on local effective field theory and 
renormalization theory, as reviewed recently in Ref. [35j | . 

A possible strategy for avoiding such complications has 
been suggested by Lifshitz and Pitaevskii |3a| : If one ac- 
cepts the Landau postulate, one may replace the true, 
singular interparticle interaction potential by a smooth, 
well-behaved one, with the requirement that the substi- 
tute shares the same value of the scattering length a. If 
then the result of the calculations, performed with the 
auxiliary potential, can be expressed entirely in terms 
of a, this result is the same as the one provided by the 
original interaction. 

In the following sections, we exploit this idea for com- 
puting the ground-state energy and the depletions of a 
binary condensate mixture by considering soft, nonsin- 
gular interaction potentials. Contrary to standard text- 
book approaches [s^ [s^l, we do not replace the inter- 
action matrix elements by constants, but keep track of 
their actual momentum dependence, so that the spuri- 
ous ultraviolet divergence is avoided altogether. Only at 
the end of the calculation, we identify the various con- 
tributions to the ground-state energy of the binary gas 
with certain terms of the Born series for the three scat- 
tering lengths involved. This allows us to cast our results 
into the desired form, enabeling a comparison with the 
known single-species results. Thus, we do neither invoke 
the pseudopotential ^ nor its refined form , but work 
directly with nonsingular interaction potentials. 

We proceed as follows: In Sec. ^ the quadratic low- 
temperature approximation to the full Hamiltonian of 
the binary gas is obtained by adapting the reasoning cus- 
tomarily applied in the single-species case. In Sec. IIIII 
this approximate Hamiltonian is then diagonalized ex- 
actly, relying on a method used by Hopfied in 1958 when 
studying the interaction of the radiation field with a po- 
larization field [s^ ]. In this way, we necessarily recover 
the already known results for the quasiparticle excitation 
spectra H, [l3 , but we also identify an energy shift as- 
sociated, for each wavevector fc, with the transformation 
to the quasiparticle form of the Hamiltonian. In order 
to obtain the ground-state energy of the system, the sum 
over all these shifts has to be evaluated; this is achieved in 
Sec. IIVI with the help of a density expansion which yields 
only finite quantities at each intermediate step of the cal- 
culation '33| . The depletions for a binary condensate are 
determined by the coefficients of the quasiparticle trans- 
formation, and calculated at the end of Sec. IIVI In Sec.lVl 
we briefly discuss our main findings. 



The Hamiltonian of a homogeneous gas consisting of a 
mixture of Ni Bose particles of species labeled by i (with 
i = 1, 2, . . .) in a cubic volume V reads 
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J2 [ dr'i'l{r)h,{r)^.^{r) 
+ \T. I I dridr2vl/|(ri)vl/t(r2) 



V JV 

xC/y(ri-r2)«',(r2)*.(ri) ,(4) 

where Uij{r) — Uji{r) denotes the real, isotropic interac- 
tion potential between particles of the type i and those of 
type j. We assume that the intraspecies potentials Uii{r) 
be repulsive; none of the potentials should admit any 
many-body bound state. The single-particle Hamilto- 
nian hi{r) is given by the kinetic energy operator for free 
particles of mass m^, hi{r) — —{h^ /2mi)Ar- The ther- 
modynamic limit will be taken at a later stage, Ni oo 
and V ^ oo such that the density = Ni/V of each 
species remains constant. The ^^(t') are the usual Bose 
field operators, obeying the commutation relations 



[vI/,(r),M'](r')] -^,,<5(r-r'). 



(5) 



Imposing periodic boundary conditions, the single- 
particle eigenfunctions of hi(r) are given by plane waves, 
y^kir) — V~^^^ exp{ikr), labeled by wavevectors k — 
^-kV^^I^ti, where n is a vector with integer components. 
The field operators then are expanded in this basis. 



(6) 



The coefficients a^fc and their adjoints ajj^. are the famil- 
iar Bose annihilation and creation operators for particles 
of species i in the single-particle state fe. As a conse- 
quence of Eq. jSjl , they satisfy the commutation relations 



(7) 



Introducing the Fourier transforms of the interaction po- 
tentials, 
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and inserting the expansion (|HJl into the Hamiltonian Q , 
one obtains 
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where the prime at the summation sign indicates restric- 
tion to those terms which comply with momentum con- 
servation, fci -|- fc2 = k[ + k'2. Note that Uij{q) carries 
the dimension of energy times volume; moreover, we have 
Utji-q) = u^{q). 

For interaction potentials with a strong short-range 
singularity, such as Lcnnard-Jones or hard-sphere poten- 
tials, the Fourier components (jHl) diverge for q ^ 0; this 
divergence can be seen as a consequence of the use of a 
single-particle basis ipkir). Therefore, keeping in mind 
the remarks made in the introduction, we perform the 
following analysis with "soft" interactions for which the 
Fourier transforms remain finite, while stipulating that 
their scattering lengths be the same as those of typical 
atomic potentials 36]. 

Under the assumption that there exists a condensate 
for each particle species, such that the single-particle 
ground state is occupied by macroscopically large num- 
bers of particles of both types, we employ the conven- 
tional quadratic approximation for the description of the 
low-energy dynamics (2^ [s^ [s^] : Firstly, the interaction 
terms of the Hamiltonian are categorized with respect 
to the number of ground-state operators Oio or ajg they 

contain. The relevant matrix elements of these opera- 

1/2 

tors should be roughly on the order of A^^ ' , and thus 

large compared to the corresponding matrix elements of 

the operators Uik or al,^ referring to excited states, and 

to sums over such elements. However, it will become 

obvious in the following that this handwaving reasoning 

requires some caution. Secondly, the "off-diagonal" oper- 
f 1/2 
ators Qi and a] g are replaced by the square roots q 

of the single-particle ground-state occupation numbers. 
This is acceptable as long as the probability distribu- 
tions for finding Nio particles in the ground state are 
broad compared to 1, so that shifts by one particle do 
not essentially change the many-particle ground state. 
Neglecting supposedly small contributions to the Hamil- 
tonian with three or more nonzero wavenumber indices, 
one finds 

H (fci - k[ )al^^ a] ^, a, ^, 
{fe} 

+ Yl H (^) -fc 



+u,j (fc) ali^a. ^/N.oNjo 
+u,j (fc) a] fcfl^ V'AWV^ 

-fw,y(0)aLa.fciV,o . 
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fe^oaifcaifc can be 



Now the relations Nio — Ni — '^k^o'^l 
employed to express the ground-state occupation num- 
bers Ni in terms of the total particle numbers Ni of each 



species. Again neglecting terms which are considered as 
small within the above categorization scheme, this leads 
to the cancelation of contributions with zero momentum 
transfer: 

C^i - fe'i ) fci a] ^, ^, 

{k} 

~ u^J{0)N,NJ + Y u,j{k)^N,Nj 
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.(11) 



Restricting ourselves from here on to only two particle 
species 1 and 2, within the quadratic approximation the 
low-energy Hamiltonian takes the form 



Hqa = Ao + J2 Hk 

k^O 



(12) 



where the superscript 2tt at the summation sign is meant 
to indicate restriction of the sum to one fc-halfspace, and 



Hk = J2[^'''{"-lk4-k + '^ik"-i-k) +(^ikalka,k 
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Dk (a{ k4-k + Oi k'^2 -fc + fc«2 fc + ^1 fc4 
(fe 
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The symbol { k —k ) denotes contributions of the 
same kind as appearing in the first two lines of this defi- 
nition, with k replaced by —k. Moreover, we have intro- 
duced the abbreviations 



Ao = ^vn{0)Ni + ^V22{0)N2 



Bik 


— Bi _fc 


Ci k 


= Ci -k 
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- D-k 



^V^2{Q)^/NiN2 



= -jV^^{k) , 



with 



h^k^ 

Irrii 

Vij{k) = ^nin.jUij{k) 



(14a) 

(14b) 

(14c) 
(14d) 

(15a) 
(15b) 



which will be used extensively in the following. Thus, 
within the quadratic approximation the Hamiltonian of 
the Bose gas mixture is given by a sum of commuting 
parts Hk, quadratic in the a-operators, which can be 
diagonalized separately; the interactions couple particles 
of both types with momenta fik and —hk. 



III. QUASIPARTICLE FORM OF THE 
HAMILTONIAN 



In order to convert Hk into a form which corresponds 
to noninteracting quasiparticles, we consider excitations 
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characterized by a wavevector fc ^ 0. Such excitations 
are described by hnear combinations of annihilation op- 
erators for particles with momentum hk and creation op- 
erators for particles with reverse momentum —hk, 



+ mka.2k + Vika2-k ^ (16) 



so that Ui k and Vi k quantify the contributions of species 1 
to the zth quasiparticle excitation, and Hik and Vik those 
of species 2. As a consequence of the symmetry of the 
coefficients (|14|) . we also have Uik = Ui-k, etc. We then 
require the commutation relations [s^ 



^ik^Hk 



(17) 



in analogy to the relation fulfilled by the Hamiltonian 
and the annihilation operator of the simple harmonic os- 
cillator. Inserting the ansatz (|16|l and the expression H13|) 
for Hk into these relations, utilizing the symmetry of the 
coefficients under interchange of k and —k, and com- 
paring cofficients of identical operators, leads us to the 
eigenvalue problem 



QkVik = Si{k)vik , 
where the 4x4 matrix Qk is given by 



Ci k 




Dk 


-Dk \ 


2-Bifc 


— Cl k 


Dk 


-Dk 


Dk 


-Dk 


C2k 


— 2i?2fc 1 


Dk 


-Dk 


2i?2fe 


— C2k 1 



and 



Vik 



{uik Vik fJ-ik Vik) 



(18) 



(19) 



(20) 



Since Qk is not hermitian, the eigenvalues £i(fc) are not 
necessarily real. For continuing the analysis, we first as- 
sume that only real eigenvalues occur, so that also the 
components of the eigenvectors Vi k can be chosen as real 
numbers, and will state the physical condition for the 
reality of the eigenvalues later (cf. Eqs. and (EZJl). 

The specific structure of Qk can be exploited with the 
help of the matrix 



/O 1 

10 

1 

VO 1 



(21) 



Since — 1, and hQkh = —Qk, multiplying the eigen- 
value problem (|18|l from the left by h leads to 



Qkhvik = -ei{k)hvik 



(22) 



Hence, if Vk = {uk Vk l^k ^k) is an eigenvector with eigen- 
value -|-e(A;), then hvk — {vk Uk Vk /^fc)* is an eigenvector 
with eigenvalue —e{k). Therefore, the eigenvalues occur 
in pairs of opposite sign: We may fix the indices such 
that 



em 

£i{k) 



-£i(fc) 
~£2{k) , 



(23) 



while stipulating that £:i(fc) and £2{k) be positive. 
Introducing the further matrix 

10 
0-100 
10 
-1 



(24) 



which obeys g = — g^^, it is easy to see that gQk is 
hermitian. Hence, two eigenvectors belonging to differ- 
ent eigenvalues of the problem (|18|1 are orthogonal with 
respect to the bilinear form mediated by g: 



hk^k 



= for i ^ j 



(25) 



where the first equality is an immediate consequence of 
the definitions (|16l) and the commutation relations {Tj) 
obeyed by the bare particle operators. In addition, some 
lines of elementary algebra reveal that the eigenvalues 
ei{k) of the matrix Qk are real — which we do assume 
here — if and only if the matrix gQk is positive semi- 
definite. Hence, we also have 



ei{k)v\^gvik = v\^gQkVik > 



(26) 



Therefore, normalizing the eigenvectors Vi k such that 
\vl^.gvik\ = 1 (note that it is not possible to change the 
sign of v\j^gVik by normalization), we thus find that the 
sign of v\f^gVik is fixed by the sign of the corresponding 
eigenvalue: 



'i k^"jk 



'"^k9VJk 

sign(ei(A;)) % . 



(27) 



Multiplying the basic commutation relation (|17|l from 
the left by bl/^, then multiplying its adjoint relation from 
the right by 6^^, and subtracting, immediately gives 



blkKk^Hk]=Q 



(28) 



again relying on the proposition that Si{k) be real. Since 
Hk depends quadratically on creation and annihilation 
operators for the bare particles, and the transforma- 
tion H16|l is linear, this operator therefore has to have 
the form 



Hk=^ftkh\fjD^,^ 



gk 



(29) 



when expressed in terms of 6-operators, where fi k and gk 
are real numbers. Once again invoking the relation H17|l 
then yields the identities 



(30) 
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which, in the hght of the commutators H27() . allows us to 
identify fn^ with the absolute value of £i{k): 



4 

El 

i=l 



gk 



(31) 



In order to exploit these algebraic deliberations for the 
construction of proper quasiparticle annihilation and cre- 
ation operators 6^^ and h\y, pertaining to the two-species 
Bose gas, we now have to make sure that these latter 
operators satisfy the Bose commutation relations 



5k,k' 



hk^b^k']-[b\kAk]=^ 



ik ' ^jfc'] 



(32) 



In view of the orthonormality-type relations (|27|l . this 
is easy to achieve: The two eigenvectors Vxk and V2k 
associated with positive eigenvalues £i{k) and e2{k) al- 
ready provide a commutator with the required magni- 
tude and sign, and therefore directly give rise, by means 
of Eq. to quasiparticle annihilation operators. 



for ; 



1,2 ; 



we also write 



£i{k) = £i{k) for i = 1, 2 



(33) 



(34) 



Next, inspecting i = j = 3 and i = j = A, the other two 
eigenvectors v^k — hvi k and v^k — hv2 k produce a neg- 
ative sign on the r.h.s. of Eq. (|27ll . Therefore, these eigen- 
vectors are associated with quasiparticle creation opera- 
tors; evidently, we have 



with positive quasiparticle energies 



|£3(fc)H£i(fe) , |?4(fc)| =e2(fe) . 



(35) 



(36) 



Note that ei(fc) = ei{—k) for i = 1, 2, as a consequence of 
the symmetry of the coefficients (|14|l . Finally, relabeling 
the operators in Eq. (|31|l according to these prescriptions, 
rearranging, and combining all emerging c-numbers into 
energy shifts a{k), allows us to cast the operators 
into the standard forms 

Hk = a{k)+ei{k)b\^b^^+e2{k)bl^b^^^ 

+ {k^ -k) . (37) 

In this way, the above construction directly leads to the 
desired oscillator Hamiltonian for noninteracting quasi- 
particles, provided we restrict ourselves to conditions 
such that the eigenvalues of the problem (|18|) are real (or, 
equivalently, those of the matrix gQk are non-negative), 
and we select the two positive eigenvalues ei(fc) and 
£2{k): Although the eigenvalues of the problem (|18|l oc- 
cur with both signs, negative quasiparticle energies are 
strictly excluded by the algebra. 



Of particular significance for our purposes are the con- 
stants a{k) occurring in the quasiparticle Hamiltoni- 
ans Hk'. For each fc, transformation of the original op- 
erator (|13|l to the normal form (|37|l necessarily involves 
such an energy shift. Since the ground state of the gas 
mixture is characterized by the absence of all quasiparti- 
cle excitations, the ground-state energy Eq is then given 
by the sum over all these shifts, plus the constant Aq 
already appearing in the quadratic approximation \12\ : 



Eq — Aq 



Y^aik), 

k^O 



(38) 



with Aq as stated in Eq. (|14a|) . In order to determine the 
individual shifts a{k), we substitute in Eq. H37|l the orig- 
inal a-operators for the 6-operators. Comparison with 
Eq. 1131) then yields 

aik)^~eiik){vjk + '^!k)-e2ik){vlk + '^L) , (39) 
together with 



Cifc = £i(fc) (Wifc 



^ik) 



2 k 



-£2{k) {i 

C2k = £i{k){filk + '^lk) +£2{k){^l, 

Finally, utilizing Eqs. H27() in the form uf,^ 
vl'u = 1, we obtain the desired relation 



^2k 



.(40) 



^4k - 



a{k) 



(ei(fe)-f £2(fe)-Cifc-C2fc) 



£2-ti - Vii - t2 



V22) 



(41) 



The expression on the right-hand side does no longer con- 
tain the transformation coefficients. Thus, for computing 
the ground-state energy according to Eq. , the knowl- 
edge of these coefficients is not required. 

Before computing the depletions of the single-particle 
ground state, we have to invert the previous ansatz 116|l . 
and to express the bare a-operators through the quasi- 
particle 6-operators. Writing 



f hk \ 



'\ -k 
hk 



a-k 



\h-kj 



a 



"ifc 

t 

1 -k 

"'2 k 

t 



(42) 



Va2-fc/ 



the result of the above construction takes the form 

bfc = VfeOfc , (43) 
with the transformation matrix 



14 



( Uik Vik Mlfe 1^1 k 

vik uik vik Mifc 

U2k V2k fJ'2k 1^2 k 

\v2k U2k V2k tJ-2k 



(44) 



The commutation relations 1)27(1 are then written as 

VkgV^ = g , (45) 
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which imphes Vi^gV^g 
V^' = gVlg 



1, or 




(46) 



This inverse transformation now allows us to state a con- 
venient expression for the depletions: Since the many- 
particle ground state |gs) of the two-species gas is the 
quasiparticle vacuum, 



6.fc|gs)-0 forfe^O 



the occupation numbers Ni k of the single-particle state 
with wavevector k for the real particles of species i, when 
the system is in its ground state, are given by 

(Nik) = {gs\a\^aj^^\gs) 

= (gs|(ui fc^I fc - t^i fc^i _fc + "2 fc^a fe ~ ^2 fe^a -fc) 
X (uifc^ife - vikbl_k + U2kh^k - 'V2khl_k)\&i) 
= ^Ik+vlk , (48a) 
^Ik + ^lk- (48b) 



{N, 



2k 



Therefore, the depletions are obtained from the sums 



N2~N20 



k^O 
k=^0 



n k 



'Ik 



'-'2k 



'2 k 



(49a) 
(49b) 



which will be evaluated in the following section. Con- 
trary to the energy shifts (|41|l , the depletions do depend 
explicitly on the transformation coefficients. 

The solution of the eigenvalue problem H18|l is elemen- 
tary: Straightforward algebra leads to the eigenvalues 



4(fc) 



-4(fc)+4(fc) (50) 
±v/(£.\(fc) - eUk)f + l6h{k)t2ik)vUk)\ , 



where here and in the following we write e+(A;) for the 
larger of the two eigenvalues ei(fc) and S2{k), and e-(fc) 
for the smaller; moreover, Ssi{k) denotes the familiar 
quasiparticle energies of a single-species Bose gas, 



4(fe) =i,(fc) [U{k) + 2vu{k)] 



(51) 



As anticipated, with each eigenvalue of the problem H18() 
also its negative appears; however, we only have to con- 
sider the positive ones. This spectrum (|5Uf) has alr eady 
been obtained through different reasoning in Refs. 11, Ea. 
In the limit of vanishing interaction between the two 
species, vi2{k) = 0, one sees that s+{k) and e_(fe) be- 
come equal to the larger and the smaller of the single- 
species eigenvalues Ssi{k) and £32(^)7 respectively. Thus, 



without interspecies interaction the single-species solu- 
tions are recovered. "Switching on" Vi2{k) makes the 
upper and lower spectral branch repel each other: If 
the single-species spectra cross for some finite fco 7^ 0, 
so that £si(fco) = £s2{l^n) = £0, one finds ej-(fco) ~ 
Eq ± 2[ti(fco)^2('i!o)wi2(^o)]^^^, so that the crossing of 
£+(fc) and £-(fe) is avoided as a consequence of the in- 
terspecies interaction. 

For fc ^ 0, both spectral branches become phononlike, 
i.e. proportional to h\k\ [sUTsf: 



4(fc) 



cWk'^ 



(47) with sound velocities c+ and c_ given by 



'11 



^22 



± 
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Here, 



Vti{0) _ Uii{0)ni 



mi 



mi 



(52) 



(53) 



(54) 



denote the squares of the well-known sound velocities for 
the single-species gases; in addition, we have introduced 



^'12(0) _ -"^2(0)^11^2 



-12 



miTO2 



mim2 



(55) 



It is remarkable that the sign of Ui2{r), and thus the 
answer to the question whether the interspecies interac- 
tion is attractive or repulsive, does not enter into the 
spectrum (|5(JII . In contrast, the intraspecies interaction 
has to be repulsive in order to guarantee the existence of 
a homogeneous single-species condensate. 

While the existence of a honogeneous binary conden- 
sate requires that both £+(fc) and £-{k) be real, the 
smaller eigenvalue S-{k) actually is real only as long as 
the interspecies interaction is sufficiently weak, that is, 
under the condition that both 



el,ik)>0, sl^ik)>0 



(56) 



and 



4vUk) < [ii(fc) + 2«n(fc)] Mk) + 2v22{k)] (57) 

be satisfied. Otherwise, our stationary approach does 
not have a physically meaningful solution. The formal 
appearance of complex eigenvalues is related to an insta- 
bility leading to phase segregation: When £_ (fc) touches 
the |fc|-axis at some /cq, the corresponding mode can 
grow without cost of energy; the slightest excitation then 
makes the system collapse into inhomogencous struc- 
tures. The condition (|56|l implies that it is not possible to 
stabilize an unstable single-species condensate by adding 
particles of another species. Since the kinetic energy con- 
tributions ti{k) grow like fc^, while the potential contri- 
butions vanish for large the stability condition H57|l 
is most likely to be violated for small momenta, or large 
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wavelength-excitations. Invoking the lowest-order Born 
approximations 



.(0) 



nii 



2Trh 



(58) 



to the s-wave scattering lengths Uij for scattering events 
between particles of species i and j, where 



niij = mim.j / [rrii -|- rrij) 



(59) 



denotes the respective reduced mass, this stability con- 
dition (|57|l reduces for fc ^ to 



mi2 



< 



(0) (0) 
'^ll °22 

miim22 



(60) 



This inequality parallels the one known from mean-field 
studies of binary condensates 0, 0| . 

The transformation coefficients determining the quasi- 
particle operators 6^ and b_ follow from 



u\ - (£±+<i)2(4-£s2)^2 4 ' 

vl = (£± - <i)2 (4 - 4) i2 , 

4 = (£± + t2f (4 - 4) ^1 A 

4 = (£±-t2)'(4-4)ii4 , 



(61) 



where we have omitted the index k for clarity, and the 
normalization factors and x_ are given by 



1 



(62) 



so that the proper commutation relations , 6j_] = 1 
are satisfied. The explicit knowledge of the coefficients 
v± k and i>± k will be required in the following section for 
evaluating the depletions according to Eqs. (|^ . 



The evaluation of the sum over these shifts a{k) is notori- 
ously difficult even in the case of a single-species conden- 
sate. If one replaces the fc-dependent interaction matrix 
elements by a constant, the sum becomes ultraviolet di- 
vergent and requires careful regularization (4Qj : in the 
context of the pseudopotential method, this regulariza- 
tion is achieved by the derivative operator appearing in 
the proper pseudopotential In standard textbook 

treatments some delicate technical details of this pro- 
cedure are tacitly swept under the carpet through the 
identification of the full scattering length with only its 
first Born approximation |3^l37l|. However, the artificial 
divergence, and the subsequent need for regularization, 
can be avoided altogether if the actual fc-dependence of 
the interaction matrix elements is kept. Here we fol- 
low an approach recently developed in Ref . |39j | , and em- 
ploy a density expansion for computing the sum in the 
low-density limit: Starting from Eq. H64|l with the ab- 
breviations H15|l. somewhat tedious but straightforward 
calculations yield 



E 

k=iO 



^a(fe) 
da{k) 



= 0, 



E 

k^O 



dfiidrij 



ni ,712—0 



ni .712—0 



7li ,712—0 



= 



4i: 



2 n'kyAm, 



(65) 
(66) 
(67) 



Now we observe that the next-to-leading term a-j^ of the 
Born series 



(0) , (1) 



(68) 



for the various s-wave scatt ering le ngths is given, in the 
limit of large volume V, by 



3rin g ie r 

mm 



IV. GROUND-STATE ENERGY AND 
DEPLETION 

From the definition (|14a|l and the Born approxima- 
tions (|58|l we infer that the first contribution to the 
ground-state energy (|38() can be written as 



(0) 

Qii 
2mii 



(0) 
«22 "2 

2to22 

JO) 



No 



(63) 



mi2 



VN1N2 



As already emphasized, when performing the transfor- 
mation from the Hamiltonian l|13() to the quasiparticle 
form (|37|l . each wavevector k furnishes the additional 
contribution H41|) . namely. 



a(k) 



e+ik) + e-{k) -Y,Mk) + vuik)] 



(64) 



(!)__/ TTiij \ ^ An ^ 4 (fc) 



(69) 



which implies 



E 



d^a{k) 



, duidrii 



2TTh^ 



(70) 



Therefore, to second order in the densities the sum over 
the energy shifts is given by 



^a(fc) ~ 2TTh^ 



(1) (1) 

Ni + -22-^N2 



2mii 



2mo 



(71) 



(1) / 

mi2 



VN1N2 



Adding Eqs. (|63|) and H71|) . we thus find the approximate 
energy density for the ground state of the binary conden- 
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sate, 

El 
V 



27rr 



2mii 



2m22 



mi2 



(72) 



At this point, a major weakness of the quadratic ap- 
proximation employed for the reduction of the orignal 
Hamihonian becomes apparent: Although the con- 
tribution (|63|l arises from the interaction term with all 
four wavevector-indices equal to the condensate index 0, 
and the contribution ifTTl) stems from terms supposed to 
be "quadratically small" , they are of the same order in 
the densities. The intuitive categorization of the interac- 
tion terms in the original Hamiltonian @ with respect to 
the number of ground-state operators they contain thus 
does not reflect a systematic ordering. The "quadratic" 
contributions to the Hamiltonian are small compared to 
the leading one, Aq, only if the Born series for the scat- 
tering lengths converge rapidly, so that 



i4JUi4°^i 



(73) 



implying that the first terms a-^'' already provide good 
approximations to the full scattering lengths fly . 

If one accepts the approximation (|72|l . one can derive a 
criterion for phase segregation by comparing the ground- 
state energy (|72|) of a homogeneous mixture of Ni and 
A^2 Bose particles in a volume V with that of two single- 
species condensates occupying separate volumes Vi and 
V2 with A^i and N2 particles, respectively, while main- 
taining a constant total volume V ^ Vi + V2- Neglecting 
surface effects, the mixture is stable as long as its ground- 
state energy is less than the lowest possible value the sum 
of two separate ground-state energies can attain for any 
partition V1+V2 — V . Performing the minimization, 
this stability requirement is translated into the condition 



4n4y V J-n+4V) (4°2^ + - 



,(1) 



TO12 



miiTO22 



(74) 



Since we require the inequalities H73|) for consistency, this 
stability condition is practically the same as the condi- 
tion deduced from the requirement of real quasipar- 
ticle energies. 

Continuing with the density expansion of the ground- 
state energy, the next-order contribution turns out to be 
of intriguing complexity. While our general scheme is 
fully capable of treating the general case, and one may 
use this scheme to generate the corresponding terms with 
the help of symbolic algebraic manipulation software, 
we restrict ourselves here to the experimentally relevant 
case of equal masses, mi = m2 = m or ti = t2 = t, 



as corresponding, e.g., to a binary condensate consist- 
ing of atoms of a single isotope in two different spin 
states 0, 0, S S I3- The resulting expressions then 
simplify considerably: The quasiparticle energies e±{k) 
acquire a form which is strongly reminiscent of the single- 
species expressions (|5T|l . 



4(fc) = t{k)[t{k) + 2u±{k)n] . 
Here we have introduced the auxiliary functions 



u±{k) 



ClUll(fc) + C2M22(fe) 



(75) 



(76) 



±J[ciUii{k) - C2U22(fc)]2 -t- AciC2ul^{k) 



which depend not only on the three interaction poten- 
tials, but also on the concentrations ci and C2 of the 
individual species. We also write n — ni -\- n2 for the to- 
tal density, implying both ni = and ci + C2 = 1. For 
vanishing interspecies interaction, Ui2(fc) = 0, one finds 
that M+(fc) and U-{k) become equal to the maximum 
and the minimum of ciuii(fc) and 02^22(^)7 respectively. 
The explicit appearance of the total density n in Eq. I|75|) 
greatly faciliates the following calculations. 

It seems to be tempting now to go on with the Taylor 
expansion (|55|l - (|67|l in powers of the density, and to 
compute 

d^a[k) i u\{k) + u^_{k) 



E 



n=0 



2 ^ 



t^ik) 



(77) 



but this expression obviously becomes infrared divergent 
in the thermodynamic limit, and therefore is of no use 
here. This finding reflects the fact that, analogous to 
the single-species case, the proper expansion parameters 
for the ground-state energy of a dilute binary Bose gas 
are not the parameters riia^j themselves, but their square 
roots. Correspondingly, the expression 

d d^a{k) 



dr 



(78) 



lim 



V 



.. *o (27r)3 
remains finite 



dfc 



iitu+f 



\5/2 



(^2 + 2tu-n) 



5/2 



(^2 + 2tu+n) 

'3^. For evaluating these integrals, we 
introduce a dimensionless momentum variable x by de- 
manding 



t = ?j.2fe2^2m : 



unx 



(79) 



with a constant u carrying the same dimension as u±{k), 
energy times volume. The limit n ^ can then be taken 
under the integral, implying u±{k) — > w±(0) there, and 
thus leading to 

a2a(fc) 



E 

{2iTn) 



d 



^{2mufl^ 3w 



(80) 



d-E 



(a;2 + 2w±(0)/m) 



5/2 
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With the help of the identity 

.,3 



dx ■ 



X 



lo (a;2 + cf/^ iVc 
one then finds 

d d^a{k) 



E 



5/2 



Introducing the scattering length-hke quantities 

.2-"±(0) 



1 
2 



(81) 



(82) 



(83) 



r n^'^'^ 4-r n''^'^ 

ciaii + C2a22 



± 



(0) 



(0) 

C2ay 



(or 



which, however, again depend exphcitly on the concen- 
trations ci and C2 of the individual species, the next con- 
tribution to the ground-state energy acquires the form 



E 



d d^a{k) 



d^fn dv? 



22 
— ? 

53 



5/2 



n=0 



128 
157^ 



E« 



(0) 



(84) 



Hence, adding this result to the previous contribu- 
tion (|7^ . we obtain for the ground-state energy density 
of a homogeneous binary condensate, consisting of two 
types of Bosons with equal masses, the expression 



27rr 



(0) 
'22 



'■22 



■ 2ciC2 (^a'{ 

E 



(0) 



'l2 



128 ■r-^ (0) 



{0)''- 



(85) 



In the limiting case of vanishing interspecies interaction 
one has 



AO) 



givmg 



V 



r (0) (o)i 
max{cia\^' , C2a22 } , 

• r (0) 



(0) , (1) , (0) 128 



(86) 



15V^ 



Since we are relying on the assumption that the Born 
series for the scattering lengths converge rapidly, so that 
a ~ a\^^ + a[l^ ~ a[^\ this formula l|57|) reduces prop- 
erly to the expected sum of two Lee- Yang-type expres- 
sions (121) ■ 



Finally, the calculation of the depletions according to 
Eqs. H49|) proceeds along similar lines. For species 1 we 
have to evaluate 



(4 



+ k 



Vlk) 



with the transformation coefficients v± k read off from 
Eqs. H61I) . Restricting ourselves again to the case of equal 
masses, and employing once more the continuous dimcn- 
sionless variable x introduced in Eq. H79|l for converting 
the sum over the discrete wavevectors k in the thermo- 
dynamic limit into an integral, we find 



'±k 



4e±t {2ei - 



2) 



(89) 



U± - C2U22 



2 2m4 



ClWll - C2U22 



■ u±/u 



+ 2u±/i 



Taking the limit n ^ in the same manner as above 
under the integral, we obtain after some manipulations 
the desired depletion 



V 



analogously, 

N2-N20 
V 



E 



(90a) 



(0) (0) 

a'^' - C2ay 



O (0) (0) (0) 

2a^ — cia\i — 02022 



3V^ 

E 



(90b) 



JO) 



(0) 

cia\{ 



2a 



(0) 



(0) 



(0) 

C2ay 



The well-known single-species expression is easily recov- 
ered in the limit of vanishing interspecies interaction, as 
it should: Utihzing Eqs. the above sums reduce to 

a single term, and we are left with 



V 



30F 



(91) 



CONCLUSION 



The preceding calculations explicitly refer to soft in- 
teraction potentials Uij{r) with a well-behaved Fourier 
transform, possessing a rapidly converging Born series 
for the respective scattering length, as is the case for 
wide and shallow potentials. In that case the first Born 
(87) 

approximations af^'' practically exhaust the full Born se- 



ries. Since the functional dependence of at least the lead- 
ing terms of the density expansion of the ground-state 
energy on the scattering lengths is the same for both 
soft interactions and hard ones with a short-range sin- 
gularity |36||] (see also Ref. for a recent discussion of 
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this point), we may free ourselves from this restriction 0^2 < 011022 for a homogeneous binary condensate of 



by replacing a[^\ or a|°^ + a\^\ by the full scattering 
lengths: Given a binary condensate of two species with 
equal masses m, total density n, and species concentra- 
tions ci and C2, interacting via some potentials Uij{r) 
with scattering lengths a^ , we define 



ClOll -I- €2022 



(92) 



±\J [ciflii - 02022]^ + 4ciC2a^2 



If the mixture is sufficiently dilute, such that 



/na^ < 1 , (93) 
the ground-state energy density is then given by 



V 



m 



Ci ail + C2 022 + 2ciC2 ai2 



+ 1^ (a+y^ + a_y^) 



, (94) 



while the depletion can be written as 



Ni~Nio 
V 



+ 



C2a22 



2a+ — cioii — C2a22 
a_ - C2a22 



2a_ 



ciaii - C2a22 



(95) 



for species 1; for the other species the depletion is ob- 
tained by interchanging the indices 1 and 2. 

The particular case — '^11022 implies a„ = 0, so 
that, according to Eq. {TSJl, one also has e-{k) = t{k), 
independent of the densities. Hence, the lower branch 
of the spectrum then is coUisionless; this phenomenon 
is of particular interest in atom interferometry |lO|. If 
one disregards the contributions to the ground-state en- 
ergy obtained in Eq. (|84|l . an argument analogous to 
the one which led to the stability criterion H74|) requires 



particles with equal masses, and thus places a system 
with a_ = right on the borderline between the regime 
of the homogeneous mixture and that of two segregated 
phases. However, inclusion of these contributions H84() 
into the analysis shifts a system with a_ = slightly 
into the regime where the mixture is stable: The contri- 
butions to the ground-state energy beyond those given 
by the mean-field approach tend to stabilize a homoge- 
neous binary condensate with a coUisionless branch of the 
quasiparticle spectrum. 

Even though derived under the restriction of equal 
masses of both species, both the result H94|l for the 
ground-state energy and, in particular, the expres- 
sion (|95|l for the depletion is surprisingly involved and 
could not have been guessed on the basis of the known 
single-species formulas. The corresponding expressions 
for the general case, which are substantially more cum- 
bersome, are obtained by following exactly the steps ex- 
plained above. 

Our reasoning clearly indicates that the elementary 
Bogoliubov theory, which works directly with the in- 
terparticle potentials, gets worse the deeper the poten- 
tials 1^ , as the missing terms of the Born series become 
increasingly more important. In the case of a single- 
species condensate, the emergence of the full scattering 
length a instead of its first Born approximation in the 
expression ^ for the ground-state energy density is jus- 
tified by means of a summation of ladder diagrams |42j | , 
or by employing effective field theory with the exact two- 
point scattering amplitude [s^ ; these devices may also be 
used in the binary case. It should be kept in mind, how- 
ever, that so far only the leading term in Eq. (|21) could be 
established with full mathematical rigor [23, whereas a 
corresponding rigorous proof for the next-to-leading term 
is still lacking. It would, therefore, be of considerable 
conceptual value to demonstrate in explicit detail how 
the missing terms of the Born series are recovered if one 
systematically takes into account those contributions to 
the full Hamiltonian which have been neglected in the 
quadratic Bogoliubov approximation, thus showing the 
consistency of the various approaches. We will present 
such a calculation in a subsequent paper |43| . 
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